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ABSTRACT 



o 

(N 

X> 

<D 

(N 

o 

u 

Oh- 

6 

(N 
> 

m 

(N 
l> 
O 

(N 
O 
O 



X 



Aims. Broad absorption line quasars (BALQSOs) are key objects for studying the structure and emission/absorption properties of 
AGN. However, despite their fundamental importance, the properties of BALQSOs remain poorly understood. To investigate the X- 
ray nature of these sources, as well as the correlations between X-ray and rest-frame UV properties, we compile a large sample of 
BALQSOs observed by XMM-Newton. 

Methods. We collect information for 88 sources from the literature and existing catalogues, creating the largest BALQSO sample 
analysed optically and in X-ray to date. We performed a full X-ray spectral analysis (using unabsorbed and both neutral and ionized 
absorption models) on a sample of 39 sources with higher X-ray spectral quality, and an approximate hardness-ratio analysis on the 
remaining sources. Using available optical spectra, we calculate the BALnicity index and investigate the dependence of this optical 
parameter on different X-ray properties. 

Results. Using the neutral absorption model, we find that 36% of our BALQSOs have N'^ < 5 x 10 21 cirr 2 , lower than the expected 
X-ray absorption for these objects. However, when we used a physically-motivated model for the X-ray absorption in BALQSOs, i.e., 
ionized absorption, ~90% of the objects are absorbed. The observed difference in ionized properties of sources with the BALnicity 
index (BI)=0 and BI>0 might be explained by different physical conditions of the outflow and/or inclination effects. The absorption 
properties also suggest that LoBALs may be physically different objects from HiBALs. In addition, we report on a correlation between 
the ionized absorption column density and BAL parameters. There is evidence (at the 98% level) that the amount of X-ray absorption 
is correlated with the strength of high-ionization UV absorption. Not previously reported, this correlation can be naturally understood 
in virtually all BALQSO models, as being driven by the total amount of gas mass flowing towards the observer. We also find a hint of 
a correlation between the BI and the ionization level detected in X-rays. 

Key words. X-rays:general - galaxies:active - quasars: absorption lines - quasars:general 



1. Introduction 

The nature of intrinsic absorption in quasars has important im- 
plications for physical models of the AGN "central engine" and, 
in general, for developing the unified model. The subclass of 
broad absorption line (BAL) quasars, characterized by strong, 
broad, and blueshifted absorption troughs in their spectra, is 
the main manifestation of the importance of outflows in the 
AGN phenomenon. While BAL objects are important in under- 
standing the properties of quasars, the nature of BALQSOs re- 
mains poorly understood and there is still no consensus about 
whether BALQSOs differ intrinsically from non-BAL quasars. 
The small percentage of BALQSOs among quasars (~ 15%) is 
generally interpreted as an orientation effect in the unified model 
(QSOs view ed at large inclination angles close to their equa- 
torial plane, IWevmann et al.lfT"991l) . However, the properties of 
some low-ionization BALQSOs appear inconsistent with sim- 
ple unified models and can be explained only by an evolution- 
ary scenario, where BALQSOs are young or recently refu eled 
quasars (e.g.. iBoroson & Meverslll992l: iBecker et alJl2000h . In 
this model, broad absorption lines appear when a nucleus blows 
off gas and dust during a dust-enshrouded quasar phase, evolv- 
ing to non-BAL quasars. Some results for the radio, such as for 
radio BALQSOs that share several properties with young radio 



sources (e.g.,|Ljuet_aI][2008) and infrared, such as the 2MASS 
Infrared survey that cont ains an unusually hi gh fraction of bright 
BALQSOs at high z (see lUrrutia et al.l l2009) seem to be in favor 
of this idea. 

Thus, the study of BALQSOs advances significantly not 
only our understanding of the structure and emission/absorption 
physics of AGN, but also helps in developing an evolutionary 
scenario of quasars. 

While the UV and optical properties of BALQSOs are more 
or less understandable and fit in with our ideas about the ex- 
pected properties of this kind of objects, the X-ray properties 
of BALs remain controversial. Theoretical model ling of BALs 
(e.g., the radiative ly driv en disk- wind paradigm by Murra y et al.l 
1995] iProga et alJ GOOO) require X-ray column densities around 
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soft X-rays generated near the central engine. Until recently, 
the number of X-ray analysed BALQSOs was rather small and, 
in many cases, the optically selected bright objects were un- 
detectable in X-rays, implying high column densities of ab- 
sorbing gas and, therefore, supporting theoretical predictions. 
The hardness ratio ana lysis of 35 X-ray detected BALQSOs by 
iGallagher et al.l d2006l) . also detects significant intrinsic absorp- 
tion (Nh ~ 10 22 - 10 24 cm -2 ). However, a sample obtained by 
cross-correlating the Sloan Digital Sky Survey and the Second 
XM M-Newto n Ser endipitous Source (2XMM) catalogues by 
iGiustini et al.l (120081) infers lower values of neutral absorption 
than found in optically selected BALQSOs samples. While this 
result is expected, given that we require that the sources be not 
just detected in X-rays, but also have enough counts to perform 
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X-ray spectral analysis, questions about the real fraction of the 
absorbed BALQSOs and properties of the absorbed gas remain 
open. Moreover, if we know that BALQSOs are objects with 
complex absorbers, outflows and ionized material (as implied 
by the definiti on and confirmed by observations of a few bright 
sources, e.g., iGruoe et al.ll200l . how correct is the application 
of the simple neutral absorbed model and how might this influ- 
ence our view of BALQSOs? 

Another important issue today is the definition of BALQSOs 
itself, which is a bit diffuse. The traditional BALs are defined 
to have CIV absorption troughs broader than 2000 kms -1 , this 
width ensuring that the absorption is from a nuclear outflow and 
effectively excluding associated absorbers. Howeve r, it could po- 
tentially exclude unusu al or intere sting BALs (s ee Becker et al.l 
2000; Hall et alj|2003) . Therefore. iTrump et al.l f2006). in their 
SDSS sample, also included in the BALQSO class BAL absorp- 
tion features at lower outflow velocities (within 3000 kms -1 ). 
As a result, a significantly higher fraction of QSOs can be clas- 
sified as BALQSOs (e.g., in the 2MASS survey, the fraction 
of BALQ SOs increases by a factor of two with the new def- 
inition, iDai et all 120081) and BALQSOs are objects with both 
weaker and much narrower absorption troughs than previously 
assu med. Despite the obvious weak points of this approach 
(see lKnigge etaDl20 08). this classification is widely used nowa- 
days. Thus, an important question arises: is there any differ- 
ence in physical properties between the "classical" and "new" 
BALQSOs? 

To explore all the questions mentioned above, we created 
the most complete sample of X-ray detected BALQSOs to date, 
based on XMM-Newton observations. 

This paper is organized as follows. We present the sample in 
Section 2. In Section 3, we describe the observations and their 
reduction; the X-ray and optical data analyses are presented in 
Section 4. Section 5 shows the results obtained from the X-ray 
spectral and hardness ratio analyses. In Section 6, we discuss our 
results and implications for different models. 

A concordance cosmology model with Hq = 
70 kms -1 Mpc~', £2a = 0.7, and Om = 0.3 is used throughout 
the paper. 



2. The sample 

Based on the material producing the BAL troughs, the 
BALQSOs are classified into high-ionization BALQSOs 
(HiBALs) and low-ionization BALQSOs (LoBALs). HiBALs 
contain strong, deep, and broad absorption troughs shortward 
of high-ionization emission lines and are typically identified by 
the presence of C IV/11549 absorption troughs. In addition to 
HiBAL features, LoBALs show broad absorption troughs in the 
low-ionization species, such as Mg 11/12798 and Al HI/11857 
lines. The rare subclass of LoBALs (~ 1%) exhibiting broad ab- 
sorption in metastable Fe II and Fe III lines is called FeLoBAL. 
About 10 to 20% of all quasars are BALQSOs. Out of these, 
about 15% are LoBALs. 

Traditionally, BA LQSOs are chara cterized by their 
BALnicity index (BI, Wev mann et al.lfl99ll) . which is a mod- 
ified velocity equivalent width of the C IV. The "traditional" 
BALQSOs have BI > Cfl The "extended" definition use s the 
absorption index (AL lHall et all 120021 ITrump et al] f2006) and 
includes BALQSOs that have absorption at outflow velocities 
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Fig. 1. Redshift distribution of the BALQSOs in our sample. 
The dotted and dashed lines correspond to the HiBALs and 
LoBALs, respectively. The upper solid line shows the distri- 
bution for the total sample of 88 sources. Most of the sources 
of both types have z around 1.9, and constitute 45% of all the 
sources. 



within 3000 kms -1 of the emission-line redshift. We used both 
kinds of objects to create our sample. 

We cross-correlate the 2XMM catalogue with the NASA 
Extragalactic Database (NEDfl the largest extragalactic sources 
database. We preselected objects with Galactic latitude greater 
than 20 degrees and high quality X-ray data corresponding to 
flag zero in the 2XMM catalogue. The queries to NED were 
launched in two different modes: either directly as a socket query 
or with a batch file. The condition for source matching was to 
fall within a distance of less than 5 arcsecs or 5 times the error in 
the position of the X-ray source queried. We finally stored only 
sources of known redshift. The detailed description of the cross- 
correlation and algorithm used is discussed in Gil-Merino et al. 

(2010, in preparation). 

In this way, we include SDSS BALQSOs by ITrump et ail 

(120061) and sources from individual random observations. In 
addition, we cross- correlate the 2XMM catalogue with other 
BALQSOs samples dShen et al]l200& iGibson et al.ll2009l) based 
on the SDSS DR5 catalogue dSchneider et al.ll2007l) . These re- 
cently published catalogues have not yet been included in NED 
and, therefore, are missed in our main cross-correlation proce- 
dure. 

Our final sample consists of 88 sources (70 HiBALs and 18 
LoBALs) and spans the redshift range 0.15 < z < 5.8 with a 
peak at z ~ 1.9 (see Fig . UJ> . Among the selected BALQSOs, we 
performed a full X-ray spectral analysis of 39 sources and an X- 
ray hardness ratio spectral analysis of the remaining 49 sources. 

The details of each source in our sample are presented in 
Table [T] The first and second columns list the NED and 2XMM 



1 That is, at least 2000 kms ' wide C IV absorption trough, 
blueshifted by > 3000 km s' 1 with respect to the C IV line. This defini- 



tion is formally detailed in Sect. 14.2 



2 The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 
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names, respectively. In the next columns, the redshift, the BAL 
type, the value of BALnicity Index (see details in Sect. 14.21) . an 
abbreviation for the performed X-ray analysis, and a number of 
the individual XMM-Newton observations are reported. The last 
column shows the notes on individual sources, i.e., details about 
the calculation of the BI or the reference from which we gathered 
information about the object. 

3. Observations and data reduction 

3.1. X-ray data 

For the full spectral analysis, we selected sources with > 70 
counts detected by the EPIC-pn instrument in the 0.2 - 10 keV 
energy band. For each observation, the X-ray data were retrieved 
from the XMM-Newton Science Archive (XSA) and reduced 
with SAS v8.0.0 and the latest calibration files. The background 
light curves at energies above 10 keV were used to filter the 
data and remove high background flaring periods. We used the 
eregionanalyse task for source region optimization and max- 
imization of the signal-to-noise ratio. The evselect tool was 
used to extract the spectrum and background region, which was 
defined as a circle around the source, after masking out nearby 
objects. We extracted the spectra of each source for the pn, 
MOS1, and MOS2 detectors using only events with pattern 0-4 
(single and double) for the pn and 0-12 for the MOS cameras. 
All spectra were extracted in the 0.2 - 10 keV band, EPIC be- 
ing the most accurately calibrated. For each spectrum, we gen- 
erated a redistribution matrix file (RMF) and ancillary response 
file (ARF) using the rmfgen and arfgen tasks, respectively. 

To improve statistics, MOS 1 and MOS2 source spectra were 
combined into a mean MOS spectrum by summing the counts 
from the channels with the nominal energy range using our own 
perl script. The background spectra and calibration files were 
merged in the same fashion. We combined the spectra by weight- 
ing them for the exposure time of the individual detectors and 
taking into account the value of the backscale parameter. If 
the source was observed more than once at compatible off-axis 
angles, we merged the pn, mean MOS spectra, and the calibra- 
tion files from one observation with the corresponding spectra 
and files from other observations. This approach allowed us to 
measure the time-averaged spectral properties (i.e., the flux and 
spectral shape), ignoring possible time variability of the sources. 

3.2. Optical data 

For sources with SDSS counterparts (69 out of 88), we retrieved 
optical spectra from the SDSS databas^l Spectra for the remain- 
ing 19 objects were either provided by authors or reconstructed, 
using a special tool, from scanned graphs of the corresponding 
articles (see details in the last column of TableQ]l. 

4. X-ray and optical data analysis 

4.1. X-ray Analysis 

A full spectral analysis was performed for a sample of 39 sources 
with the highest quality statistics, leading to proper values of 
r and Nh in the source rest-frame. We used XSPEC vl2.4.0 
dArnaudll9 96) to perform the spectral analysis. All spectra were 
binned to a minimum of 3 counts/bin. We performed this min- 
imal grouping to avoid channels with no counts, i.e., the spec- 
trum is essentially unbinned and no spectral information is lost. 

3 http://cas.sdss.org/astro/en/tools/search/SQS.asp 
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Fig. 2. . The power-law photon index F versus rest-frame 
for the analysed sample of 39 sources suitable for full spectral 
analysis. Solid black squares and red circles refer to LoBALs and 
HiBALs, respectively. Open red circles indicate HiBALs with 
BI=0. The parameter T for our objects ranges from 1.18 to 2.72, 
the majority of the sources clustering around 1.9. 



Because of the small number of counts in our sp ectra, we de - 
cided to use the Cash-statistic (XSPEC C-stat, ICashlll979l) . 
The new version of XSPEC allows C-statistic fits to data for 
which background spectra are considered. Comparison of the 
fitting results from Cash and x 2 statistics (for our best quality 
spectra) reveals that the difference is negligible. Since C-stat 
does not indicate at all of quality of the fit, we computed its 
goodness using Monte Carlo probability calculations. For this 
purpose, we used the goodness command in XSPEC, which 
performs Monte Carlo simulations of 100 model spectra using 
the best-fit model and infers the percentage of simulated spec- 
tra that had a fit statistic less than that obtained from the fit to 
the real data. The 'goodness of fit' should be around 50%, if the 
observed data were produced by the fitted model. The obtained 
values were, indeed, around this value for all our spectra. 

EPIC-pn and MOS data were fitted simultaneously using ap- 
propriate models. We, initially, fit the spectra with a model con- 
sisting of a power law with an intrinsic absorption component 
at the source redshift zphabs, and an additional photoelectric 
absorption component phabs that was fi xed at the Galactic col- 
umn density (Dickey & Lockmanlll990l) . From our model fits, 
we computed the slope of a power-law spectrum (photon index 
F), the intrinsic rest-frame neutral column density , and the X- 
ray flux and luminosity in the 0.5 - 2 and 2-10 keV bands. For 
absorbed sources, the rest-frame luminosities were corrected for 
absorption by setting to zero the Nh values in the XSPEC best-fit 
model. Spectral analysis results are reported in Table [2] 

Fourteen of the 39 BALQSOs (-36%) have low intrinsic 
neutral ab sorption, N R < 5 x 1 21 cm 2 , confirming the result 
(-36%) of iGiustini etafl (120081) . 

As noted in the Introduction, this result is difficult to under- 
stand in the context of any BALQSO scenario. For this reason, 
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we decided to study the possibility that ionized (as opposed to 
neutral) gas is significantly present in our BALQSOs, and then 
characterise the properties of this ionized absorber. If BALQSOs 
are special sources with strong outflowing ionized winds and 
ionized or partially covering intrinsic absorption, then we need 
to adopt a model that takes account of the presence of these fea- 
tures. To test this, we replaced the neutral absorber in our spec- 
tral model with an ionize d intrinsic absorbe r, modelled using the 
XSPEC model absori dDone et al.lll992l) . Owing to the low 
signal-to-noise ratio of the data, we forced the absori photon 
index to be fixed and equal to the continuum power-law index 
corresponding to F obtained from the previous absorbed power- 
law model. The temperature of the absorbing material was fixed 
at T = 3.0 x 10 4 K and its redshift was fixed to be the redshift of 
the source. 

To constrain the output parameters of the ionized hydrogen 
column density and ionization parameter ^Q, we varied the 
values of the parameter £ from 10 to 1000 erg cm s _1 in steps 
of 50 erg cm s~'0 to find best-fit values of N' H , and then we 
used these values of to constrain Each of these steps was 
checked against the 'goodness of fit' by the simulation of 100 
model spectra. Using this approach, we found border values of 
two parameters for each of our sources. Almost as expected, this 
model produces higher values of absorption for all of our sources 
(only 3 sources out of 39 still have < 5 x 10 21 cm" 2 ). This is 
consistent with a scenario un which there is a significant amount 
of gas along our line of sight to the nuclei of BALQSOs, but this 
gas is likely to be ionized in most cases. The estimated parame- 
ters for the absori model are listed in Table[3] 

We performed a hardness ratio (HR) spectral analysis of our 
lower spectral quality sources (49 out of 88). The hardness ratio 
was calculated using the standard formula HR = (H-S )/(H+S ), 
where H and S correspond to counts in the hard (2.0 - 10.0 keV) 
and soft (0.5 - 2.0 keV) energy bands, respectively. Using 
XSPEC, we calculated the HRs for a grid of power-law mod- 
els with the canonical value of F = 1.8 and neutral absorption 
in the range of 10 21 - 10 24 cirT 2 at the source redshift. Both the 
Galactic and intrinsic absorption were i ncluded in the model, de - 
tails of the calculation being present in lGil-Merino et af . (2009). 
We then compared the observed hardness ratio of each source 
with the modelled ones, and derived the amount of the intrin- 
sic absorption. The absorbing column densities and their lcr er- 
rors (propagated from count-rate errors) are reported in Table [4] 
In Fig. [3] we show the absorption distributions for our HiBAL 
(top panel) and LoBAL (bottom panel) subsamples. For each 
subsample, we plot the absorption values obtained from the HR 
analysis of our lower spectral quality sources and the fitted re- 
sults for the 39 spectrally analysed sources. 

4.2. Optical data analysis 

Following the definition given by IWevmann et al.l (119911) . we 
calculated the BALnicity index of C IV to be: 
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4 £ = L/nr 2 , where L is the ionizing luminosity of the source, n is 
the number density of the absorber, and r is the distance between the 
absorber and the ionizing source 

3 We choose this range of £ values in order to reproduce possible 
ionization states of the material, i.e., from almost neutral (£ < 50 erg 
cm s -1 ) to highly ionized (£ ~ 1000 erg cm s -1 ) absorption 
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Fig. 3. Top panel: Distribution of neutral absorbing column den- 
sities for 70 HiBALs (top panel) and 18 LoBALs (bottom panel) 
in our sample. In both graphics, the red solid lines refer to val- 
ues obtained from the full spectral analysis of our high quality 
spectra (30 HiBALs and 9 LoBALs, respectively), while black 
dashed lines refer to values obtained from the HR spectral anal- 
ysis of the remaining sources. The insert in the top panel shows 
the observed distribution of absorption for the total sample of 
70 HiBALs divided by values of the BALnicity Index. The thin 
black line refers to the 49 sources with BI>0 and thick blue line 
refers to the 21 sources with BI=0. 



where f(v) is the normalized flux (calculated from observed and 
fitted fluxes) as a function of velocity in km s , and C = 1 
at trough velocities greater than 2000 kms~' from the start 
of a contiguous trough, and C = elsewhere. In the case of 
the calculation of BI from Mg II, which is generally narrower 
than the high-ionization lines, the integral in Eq. (1) starts from 
v = kms -1 and C = 1 at trough velocities more than 1000 
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5. Results 
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Fig. 4. Top panel: The ionization level of absorbing material 
f versus ionized column densities. Solid black squares and red 
circles refer to LoBALs and HiBALs, respectively. Open red cir- 
cles indicate HiBALs with BI=0. Bottom panel: BALnicity in- 
dex versus Solid black squares and red circles refer to LoBALs 
and HiBALs, respectively. The black stars refer to LoBALs with 
BI calculated from Mg II. 



km s 1 . The BALnicity index can take on any value in the range 
kms '< BI <20,000 kms 1 . 

Because of the redshift range of our LoBALs, only for 3 
sources do we have BI measurements from C IV, and for the 
remaining objects (at z < 1.8) the Mg II BI is presented (see 
TableQ}. 



To investigate the general X-ray properties of our sample of 39 
BALQSOs (9 LoBALs and 30 HiBALs) using the highest quality 
statistics and compare the results with previous findings, we first 
measure A 7 ^ and F by adopting an absorbed power-law model 
with neutral absorption. For the remaining 49 sources, upper and 
lower limits to N'^ values can be placed using the hardness ratio 
spectral analysis. 

We do not find any correlation between the photon index 
and intrinsic absorption (Fig. [2]). The mean photon index is 
(T) = 1.87, typical of values found in X-ray analyses of radio- 
quiet type-1 AGN. We detect absorption > 5xl0 22 cirT 2 
in 25 objects (~64%). If we consider the whole sample of 88 
BALQSOs, the fraction of absorbed sources increases to 67%, 
still meaning that about one third of the sample is X-ray unab- 
sorbed. We do not find any strong difference between the neu- 
tral absorption properties of the "classical" (BI>0) and those 
of "new" (BI=0) BALQSOs (Fig.©. The Kolmogorov-Smirnov 
(K-S) test infers the maximum value of the absolute difference 
between the two cumulative distribution functions, D = 0.25 
with a probability, p, that these objects are drawn from the same 
population of 0.88. 

Taking into account the difference in spectral and HR anal- 
ysis strategies (e.g., use of^ 2 instead of C-stat, assuming dif- 
ferent energy ranges for the calculation of HR), o ur results are 
in very good agreement with lGiustini e t al. (200^) for our over- 
lapping 54 BALQSOs. 

It is known that the complex structure of a spectrum (pres- 
ence of ionized/partially covering absorption or additional fea- 
tures in the underlying continuum, such as the soft-excess or 
scattered component) can conspire to produce apparently lower 
values of the intrinsic absorption. Taking into account that 
BALQSOs are sources with complex ionized structures by defi- 
nition, we fit our 39 objects with the absori model. 

This fit of a power law absorbed by ionized material leads 
to a slight improvement in the 'goodness' with respect to the 
neutral absorbed power-law model, which is however not sta- 
tistically significant in many cases, mainly because of the poor 
quality of the data) and provides an acceptable parameterization 
of the spectra. It should be emphasized that the need for an ion- 
ized absorber (instead of a neutral one) is certainly not required 
by the X-ray data, but by the likely physical condition of the ab- 
sorbing gas. forthermore, our purpose is to study and constrain 
the possible ionized properties of BALQSOs, rather than com- 
pare the goodness of fit for different possible models. 

In general, the absori model provides higher values of 
A^j for almost all our sources. We detected absorption > 
5xl0 21 cirT 2 in 36 sources (~92%) and 7 of them showed an 
absorption ~10 23 cirT 2 . A plot of ionized column density ver- 
sus ionization parameter £ is given in Fig. [4] (top panel). There 
is a clear correlation between £ and N' H for our HiBALs con- 
firmed by both the Kendall's tau (r = 0.55 with a probabil- 
ity of p = 1.99xl0~ 5 ) and Spearman's rho (p = 0.69 with 
p = 2.41xl0~ 5 ) tests. 

Some of the LoBALs appear as clear outliers in this graphic 
and, in general, these sources have lower values of than 
HiBALs. The LoBALs also do not exhibit the stronger neutral 
absorption that HiBALs do: only in 2 LoBALs (out of 18) do we 
detect A^ ~10 23 crrT 2 , while most of the objects have lower val- 
ues (<10 22 cm" 2 ). The K-S test confirms that the LoBALs have 
a different absorption distribution than HiBALs (p = 0.067 and 
p = 5.89xl0~ 3 for the ionized and neutral absorptions, respec- 
tively). The LoBALs-outliers show the highest ionization level 
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of absorbing material (although it is poorly constrained) and the 
lowest values of absorption among the analysed sources, sug- 
gesting that the intrinsic nature of these objects must differ from 
the remaining BALQSOs. This result agrees with the idea that 
LoBALs (or at least some of them) may physically differ from 
HiBALs, and hence require the application of different models 
for describing their intrinsic properties (e.g., partial covering in- 
trinsic absorption). 

While we do not find any strong difference in the neutral 
absorption properties for our sample of HiBALs as a function 
of the BALnicity index, we do observe a separation in the case 
of ionized absorption. For our 6 HiBALs with BI=0, we observe 
mean values of (£) =110 erg cm s" 1 and (N^) = 2.9xl0 22 cm" 2 , 
while the 24 "classical" BALQSOs show (0 = 250 erg cm s" 1 
and <A^> = 8.4xl0 22 cm" 2 . The K-S test indicates only a 10% 
chance that these two samples of sources are drawn from the 
same population. 

There is a marginal tendency to measure higher ^ values with 
increasing BI (Fig. |U bottom panel). However, both Kendall's 
tau and Spearman's rho tests give weak support for a correla- 
tion. The Kendall's tau method gives r = 0.14 with an associated 
probability of p = 0.41, while the Spearman's rank correlation 
test gives p = 0.31 with p = 0.17. We also checked the cor- 
relations between BAL properties and the neutral/ionized X-ray 
absorption column density (measured through the X-ray spec- 
tral fit or HR analysis). While we do not see any strong depen- 
dence between these parameters in the case of neutral absorption 
(Fig- 13 top panel), there is an apparent trend of increasing ion- 
ized with BI for HiBALs (Fig. [5] bottom panel). Applying 
the Kendall's tau and Spearman's rho tests, we obtain r = 0.32 
and p = 0.43, corresponding to a probability that the ionized 
N' H and BI are uncorrelated of p T = 0.02 and p p = 0.018, re- 
spectively. Although the significance of the trend is only ~98%, 
it might have a very simple physical origin, for which the total 
amount of outflowing gas largely dictates the BI and the ionized 
absorption column density. 

6. Discussion 

As mentioned in the Introduction, an outstanding question about 
BALQSOs is whether they are all heavily absorbed or only 
a fraction of them are. Quite c ontroversial results were pre- 
sented by iGallagher et al ] (120061) in their Chandra analysis of 
35 Large Bright Quasa r Survey BALQSOs (all sources have 
iVg > 10 22 cm" 2 ) and iGiustini et alj d2008h in their XMM- 
Newton analysis of 54 SDSS BALQSOs, nearly half of the 
sources having < 10 22 cm" 2 . 

Our final sample spans a wide range of intrinsic absorption 
column densities (derived from the neutral absorption model), 
although the main result, that -68% of 88 B ALQSOs have 
> 5 x 10 21 cm" 2 , is in agreement with IGiustini et al] 
(20081). Thus, the BALQSOs class remains a class of, in general, 
absorbed sources. For comparison, a typic al fraction of broa d 
line AGN with neutral absorption is -3% (Mat eos et al.| [2009). 
although BALQSOs seem to be less absorbed than previously 
assumed. 

As mentioned by Giustini et al. (2008), some factors can in- 
fluence this result. One possible explanations may be the dif- 
ferent energy range used for the spectral analysis in Chandra 
and XMM-Newton data (0.5-8 keV and 0.2-10 keV, respec- 
tively). In addition, our work is biased toward the X-ray bright- 
est sources, because we searched for all known BALQSOs with 
an X-ray detection, while previous works were mainly based on 
purely optically selected BALQSOs. 




10 20 10 21 10" 1 23 10 24 

neutral N H (cm -2 ) 




0.5 - 

i i i i i i i I i i i i iii 

1 22 1 o 23 

ionized N H (cm -2 ) 



Fig. 5. Dependences of BALnicity index on neutral (top panel) 
and ionized (bottom panel) absorptions. In both panels, the solid 
black squares and red circles refer to LoBALs and HiBALs, re- 
spectively. The black stars refer to LoBALs with BI calculated 
from Mg II. Open red pentagons and black triangles refer to 
HiBALs and LoBALs, whose values of Nh were obtained from 
HR analysis. 



Does this result mean that we need to revise our ideas about 
the nature of BALQSO? The answer is not as obvious as it 
seems when we use physically more motivated approaches for 
the modelling of the inner absorption in these objects. Knowing 
the complicated nature of these sources, we used the absori 
model to constrain the ionized properties of BALQSOs. Using 
this model, we detected absorption > 5xl0 21 cm" 2 in >90% 
of them. The average value of the absorption is in good agree- 
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ment with the predictions of a radiatively driven accretion disk 
wind model bv lMurrav et alJ (1 19951) . The observed properties of 
LoBALs confirm the idea that these objects may physically dif- 
fer from HiBALs. 

Using our sample we also have the possibility of comparing 
properties of "classical" (BI>0) and "new" (BI=0) BALQSOs. 
We found that the sources with BI=0 lack significant X-ray ab- 
sorption, and, in general, this subsample hosts the lowest of 
the entire sample. In our study, the BI>0 subsample appears to 
be more X-ray absorbed and ionized than the BI=0 one, which 
is clearly consistent with any picture where the mass of out- 
flowing gas largely determines both the optical/UV absorption 
troughs and the X-ray absorption. In other words, even if we can- 
not reject the possibility that there are unabsorbed, X-ray bright 
BALQSOs, it seems that the "classical" BALQSOs (i.e., with 
BI>0) are always X-ray absorbed. 

Our results suggest that "new" BALQSOs are, likely, a dif- 
ferent class of absorbed quasars from the "classical " BAL QSOs. 
A similar conclusion is reached by Knig ge et al.l (120081) . who 
studied the QSO optical spectra in the SDSS DR3 catalogue. The 
differences in the observational properties of these two classes of 
objects may be related to the large angle between the line of sight 
and the axi s of the accre tion disk plane. In the QSO unification 
scheme of lElvisI d2000l) . intrinsic absorption lines of different 
widths form in the same disc wind and are just representations of 
different lines of sight through the outflowing wind to the quasar 
continuum source. The different properties of quasars may also 
be due to real physical differences between the "new" and "clas- 
sical" BALQSOs. In particular, the less ionized and weaker out- 
flows of "new" BALQSOs may represent the late evolutionary 
stages (i.e., dissipatio n over time) of "classic al" massive BAL 
winds (see review bv lHamann & Sabrall2004 . In our opinion, 
both factors (orientation and evolution of the flow) are likely to 
contribute to the observed differences between these two classes. 

A purpose of this paper was also to study the relationship 
between the BAL properties and X-ray absorption. For the first 
time, we find a correlation between the ionized absorption and 
the BALnicity index (i.e., the UV absorption line properties). 
Although the correlation is significant at the 98% level, our re- 
sults demonstrate that the amount of X-ray absorption and the 
strength of UV absorption are related when an ionized absorp- 
tion model is used to model the X-ray data. The physics of this 
could be as simple as the mass of outflowing gas. 

This result may help us to understand the difference between 
our neutral absorbing column density values and those of the 
Gallagher et al.l (120061) sample. If we assume a direct depen- 
dence between the BALnicity index and absorption, then in the 
samples with higher BI we must, in general, observe higher X- 
ray column densities. Although the mean value of the BI in the 
LBQS BALQSOs sample is (BI) =3437 kms" 1 , in our sample 
we have just a few objects with BI> 3500 kms _1 ((BI) =1077 
kms _I ). The neutral column densities of almost all our high BI 
objec ts are >10 22 cm" 2 and , thus, compatible with the values 
from lGallagher et al l d2006l) . 

Summarizing our results, we can confirm the idea that the 
presence of the X-ray ionized absorption is important for launch- 
ing BAL winds. While the influence of other factors on the 
properties of the wind cannot be rejected, it seems that this de- 
pendence is clear: the higher the X-ray absorption, the higher 
the UV absorption and, hence, the more powerful the outflow. 
However, the required absorption values for outflow launching 
mechanisms are not as high as has been previously supposed, 
which might actually be a selection effect of our X-ray selection 
bias. 
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Table 1. 2XMM BALQSO sample 
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Name (1) 


2XMM name' 2 ' 




BAL Type (4) 


BP' 


Analysis* 6 ' 


N. obs. (,) 


Notes (!i) 


2QZ J003814.3-273130 


2XMM J003814. 1-273 131 


2.405 


H 


597 


HR 


1 


Cr 


SDSS J004206.18-091255.7 


2XMM J004206.1-091255 


1.778 


nHi 





HR 


1 




SDSS J004338. 10+004615.9 


2XMM J004338.0+004616 


1.574 


Lo 


185 


A 


1 


Mg II 


SDSS J004341. 24+005253.3 


2XMM J004341. 3+005253 


0.834 


Lo 


91 


A 


1 


Mg II 


APMUKS(B J)B00475 1 .94-523 127.9 


2XMM J005007.3-521506 


2.422 


Hi 


1300 


A 




C 


XMMU J0103 16.7-065 137 


2XMMJ0103 16.2-065 135 


1.914 


Lo 


750 


HR 


1 


Br 


XMMU J010328.7-064633 


2XMM J010328.7-064634 


1.820 


Hi 


210 


A 


1 


Br 


SDSS J010941.97+132843.8 


2XMM JO 10941. 8+ 132844 


1.226 


Lo 


20 


HR 


1 


Mg II 


[WFM91] 0226-1024 


2XMM J022839. 1-101111 


2.256 


Hi 


7344 


HR 


1 


W 


SDSS J023 148.80-073906.3 


2XMM J023 148.8-073905 


2.512 


Hi 


50 


A 


1 




SDSS J023224.87-071910.5 


2XMM J023224.7-071909 


1.597 


Hi 


280 


A 


1 




SDSS J024230.65-000029.6 


2XMM J024230.6-000030 


2.505 


H 


2286 


HR 


1 




SDSS J024304.68+000005.4 


2XMM J024304.6+000005 


1.995 


Hi 


290 


A 






SDSS J030222.08+000631.0 


2XMM J030222.0+000630 


3.315 


H 


100 


A 


1 




SDSS J072945.33+370031.9 


2XMM J072945.3+370032 


0.969 


Lo 


37 


HR 


1 


Mgll 


SDSS J073405.24+3203 15.2 


2XMM J073405.2+320315 


2.082 


Hi 





HR 


1 


IRAS J07598+6508 


2XMM J080430.4+64595 1 


0.149 


Lo 


3892 


A 


1 


H 


SDSS J083633.54+553245.0 


2XMM J083633.4+553245 


1.614 


Hi 





HR 


1 




CXOSEXSI J084840.5+445732 


2XMM J084840.5+445733 


3.093 


Hi 


2136 


A 




St 


SDSS J08555 1.24+375752.2 


2XMM J08555 1 . 1 +375752 


1.930 


Hi 


1296 


A 


1 




SDSS J091127. 61+055054.1 


2XMM J091 127.5+055054 


2.763 


H 


2358 


A 


1 




SDSS J091914.23+303019.0 


2XMM J091914.2+303018 


1.387 


Lo 


38 


HR 


1 


Mgll 


SDSS J092104.36+302030.3 


2XMM J092104.3+302031 


3.350 


H 


40 


A 






SDSS J092138.45+301546.9 


2XMM J092138.4+301546 


1.589 


Hi 





HR 


1 




SDSS J092142. 57+513149.4 


2XMMJ092142.5+513148 


1.835 


nHi 


67 


HR 


1 




SDSS J092238.43+512121.2 


2XMM J092238.3+5 12120 


1.751 


Lo 


180 


A 


1 


Mgll 


SDSS J092345. 19+512710.0 


2XMMJ092345.1+512711 


2.164 


Hi 


1696 


HR 


1 


SDSS J093918.07+355615.0 


2XMM J093918. 1+355612 


2.047 


Hi 


80 


HR 


1 




SDSS J094309.56+481 140.5 


2XMM J094309.8+481 142 


1.809 


Hi 


65 


HR 


1 




SDSS J095548. 13+410955.3 


2XMM J095548. 1+410955 


2.308 


nH 


10 


HR 


1 




SDSS J095942.08+024103.1 


2XMM J095942.0+024104 


1.795 


Lo 


595 


HR 


1 


P 


SDSS J095944.47+05 1158.3 


2XMM J095944.4+05 1 157 


1.596 


Hi 





HR 


1 




SDSS J 10003 8. 01 +020822. 3 


2XMM J100037.9+020822 


1.825 


Lo 


267 


HR 


1 


P 


SDSS J1001 16.78+014053.5 


2XMM J1001 16.7+014053 


2.055 


Hi 





A 


2 




SDSS J100120.84+555349.5 


2XMM J100120.7+555351 


1.413 


Lo 


27 


A 


2 


Mgll 


SDSS J100145. 15+022456.9 


2XMM J100145.2+022456 


2.032 


nHi 


10 


HR 


2 


PG 1004+130 


2XMM J00726.0+ 124856 


0.241 


Hi 


850 


A 


1 


Wi 


SDSS J100728.69+534326.7 


2XMM J100728.8+534327 


1.768 


nHi 


3 


A 


1 




SDSS J101 144.33+554103.1 


2XMM J101 144.4+554103 


2.811 


H 





A 


1 




SDSS J101614.26+520915.7 


2XMM J101614.2+520916 


2.455 


H 


2401 


HR 


1 


Gr 


SDSS J101954. 54+082515.0 


2XMM J101954.6+082515 


3.010 


nH 





HR 


1 




SDSS 102117.74+131545.9 


2XMMJ1021 17.7+131546 


1.565 


Lo 


35 


A 


1 


Mgll 


SDSS J104433.04-012502.2 


2XMM J104433.0-012500 


5.800 


Hi 


900 


HR 


1 


G 


ISCLLHDS J105154+572408 


2XMM J105 154.5+572407 


2.400 


Hi 


900 


A 


12 


LH 


SDSS J105201. 35+441419. 8 


2XMMJ105201. 3+441417 


1.791 


Hi 


476 


A 


1 




RX J105207.7+573842 


2XMM J105207.4+573838 


2.730 


Hi 


4300 


A 


12 


LH 


SDSS J110637. 16+522233.4 


2XMM Jl 10637.0+522233 


1.840 


Hi 





HR 


1 




SDSS Jl 10853.98+522337.9 


2XMM Jl 10853.5+522341 


1.664 


Hi 


50 


HR 


1 




[HB89] 1115+080 


2XMM Jl 1 1816.9+074558 


1.735 


Hi 


9 


A 






SDSS Jl 11859.56+075606.5 


2XMM Jl 1 1859.7+075602 


1.760 


Hi 


3 


HR 


1 




SDSS Jl 12020.96+432545.1 


2XMM Jl 12020.9+432545 


3.548 


H 





A 


1 




SDSS Jl 12300.25+052451.0 


2XMM Jl 12300.3+052448 


3.700 


nH 





HR 


1 




UM425 


2XMM Jl 12320.7+013748 


1.462 


Hi 


415 


HR 


1 


Hu 


SDSS J112432.14+385104.3 


2XMMJ1 12432.0+385 104 


3.530 


H 


270 


HR 


1 




SDSS J113537. 67+491323. 2 


2XMMJ1 13537.6+491322 


1.982 


nH 


205 


HR 


1 




SDSS Jl 141 1 1.61-014306.6 


2XMM Jl 141 1 1.5-014305 


1.266 


Lo 


500 


HR 


1 


Mgll 


SDSS Jl 14636.88+472313.3 


2XMM Jl 14636.9+472313 


1.895 


Lo 


24 


A 


1 




SDSS J120522.18+443140.4 


2XMM J120522. 1 +443 141 


1.921 


Hi 


965 


A 






SDSS J122307.52+103448.2 


2XMM J122307.4+ 103448 


2.742 


H 





A 






SDSS J122637.02+013015.9 


2XMM J122636.9+013016 


1.552 


Hi 


787 


HR 






SDSS J122708.29+012638.4 


2XMM J122708.2+012638 


1.954 


Hi 


730 


HR 






SDSS J123637.45+615814.4 


2XMM J123637.7+615813 


2.520 


H 


447 


HR 






LBQS 1235+1807B 


2XMM J123820.3+ 175038 


0.449 


Lo 


5086 


HR 




Mg II, F 


SDSS J124520.72-002128.2 


2XMM J124520.6-002127 


2.354 


H 





A 






SDSS J124559. 59+570053.1 


2XMM J124559.6+570052 


1.656 


Hi 





HR 
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Name* 11 


2XMM name (i) 


z (j) 


BAL Type (4) 


BI W 


Analysis' 6 ' 


N. obs. (/) 


Notes (s) 


[HB89] 1246-057 


2XMM J124913. 8-055918 


2.236 


Hi 


4590 


A 


2 


J 


SDSS J132330.44+545955.6 


2XMM J 132330.7+545954 


2.208 


Hi 





A 


1 




SDSS J132401.53+032020.6 


2XMM J 1 3240 1 .4+0320 1 9 


0.926 


nLo 


217 


HR 


1 


Mg II 


SDSS J132827.07+581836.9 


2XMM J132827.3+581839 


3.139 


H 


198 


HR 


1 




SDSS J133553.61+514744.1 


2XMM J133553.7+514744 


1.830 


Lo 


340 


A 


1 




SDSS J133639.40+514605.2 


2XMM J 133639. 1+5 14608 


2.229 


Hi 


1410 


HR 


1 




SDSS J134059.24-001944.9 


2XMM J134059. 1-001945 


1.857 


nHi 





HR 


1 




HELLAS2XMM J140049.1+025850 


2XMM J140048.7+025852 


1.822 


Hi 


569 


HR 


2 


C 


SDSS J142539.38+375736.7 


2XMM J142539.3+375736 


1.897 


Hi 


104 


A 


1 




SDSS J142555.22+373900.7 


2XMM J142555.2+373900 


2.691 


H 


53 


HR 


1 




SDSS J142652.94+375359.9 


2XMM J142652.8+375401 


1.812 


Hi 


95 


HR 


1 




SDSS J143513. 89+484149.3 


2XMM J143513.9+484149 


1.886 


Hi 


20 


A 


2 




SDSS J144625.48+025548.6 


2XMM J144625.6+025549 


1.866 


Hi 


256 


HR 


1 




SDSS J144727.49+403206.3 


2XMM J144727.4+403206 


1.335 


Lo 


50 


A 


1 


Mg II 


SDSS J150858. 15+565226.5 


2XMM J150858.2+565227 


1.797 


Hi 





HR 


1 




SDSS J151729.70+001652.6 


2XMM J151729. 6+001651 


1.887 


nHi 


80 


HR 


1 




SDSS J152553. 89+513649.1 


2XMM J152553. 8+513649 


2.883 


Hi 


3144 


A 


3 




SDSS J153322.80+324351.4 


2XMM J153322.8+324351 


1.899 


Hi 





HR 


1 




SDSS J154359.44+535903.2 


2XMM J 154359.4+535902 


2.370 


Hi 


82 


A 


2 




SDSS J203941. 04-010201. 6 


2XMM J203941.2-010202 


2.065 


Hi 





HR 


1 




LBQS 2111-4335 


2XMMJ21 1507.0-432310 


1.708 


Hi 


7249 


A 


1 


Mor 


SDSS J213023.61 + 122252.2 


2XMM J2 1 3023 .4+ 1 2225 1 


3.263 


nH 


12 


A 


1 




SDSS J231850.79+002552.6 


2XMM J23 1850.6+002554 


1.591 


Hi 





HR 


1 





Notes: Col.(l): Source name; Col. (2): 2XMM source name; Col.(3): Redshift, taken from SDSS or NED; Col. (4): The BALQSO 
subclassification: Hi - HiBAL, Lo - LoBAL, H - HiBAL in which the Mgll region is not within the spectral coverage (i.e., impossible to check the 
presence of low-ionization absorption troughs), n - relatively narrow trough; Col. (5): BI - BALnicity Index (in units of km s~'); Col. (6): X-ray 
analysis used in this work: A - spectral analysis, HR - hardness ratio analysis; Col. (7): Number of individual XMM-Newton observations; 
Col. (8): Notes on individual objects : "Mg II" denotes L oBAL, which BALnicity index were compute d from Mg II absorp tion troughs; "Cr" - the 
optical spectrum was retrieved from Croom et al. (2004); "C" - optical spectra courtesly provided bv lCocchia et al. (2007); "B r" - op tical spectra 
from AXIS survey courtesly provided by Barco nTet al.1 i206"2h : "W" - BALnicity index value was taken fr om Wevma nn et al.1 dl99lh : "H" - 
reconstructe d spectrum fromlHines & W ills (1995); "St" - reconstructed spectrum from Stern et al. ( 20Q2); "P" - optical spectra courtesly 
provided by Prescott et al. (2006); "Wi" - BALnicity index value was t aken from Wills et al. dl999h : "Gr" - BALnicity index value was taken 
fromlGregg et alj (120001) : "G" - BALnicity index value was taken from Goodrich et al. ( 200lb ; "LH " - spectra were taken and reconstructed from 
i Lehmann et all d200lh : Hu" - BALnicity index valu e was taken from| Huts emekers & Lamvl d2000t> ; "F" - reconstruct ed spectrum from 
Foltz et all dl987t) : "J" - reconstructed spectrum from Junkkarin en et al.1 dl987t) : "Mor" - reconstructed spectrum from Morris et alHl99l|). The 
spectr a for rest of the objects w ere retrieved from the SDSS archive and correspond to the following catalogues: [Trump et al.1 J20 06): Shen et al] 
(2008) and lGibson et al.U2009l) . 
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Table 2. Best fit parameters from the absorbed power-law model for a sample of 39 BALQSOs with the best statistics. 
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Notes: Col.(l): Source name; Col. (2): Galactic neutral hydrogen column density in units of 10 20 cnT 2 , taken from Dic key & Lockmanl dl990h : 
Col. (3): Photon Index; Col. (4): Neutral hydrogen column density, in units of 10 22 cm" 2 ; Col. (5) and Col. (6): Logarithms of 0.5 - 2 and 
2-10 keV observed fluxes, in units of erg s~' cm 4 ; Col. (7) and Col. (8): Logarithms of 0.5 - 2 and 2-10 keV rest-frame luminosities, corrected 
for intrinsic absorption, in units of erg s _1 . 



A. Streblyanska et al.: BALQSOs 
Table 3. absori model fits to the X-ray spectra for a sample of 39 BALQSOs with the best statistics. 
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Notes: Col.(l): Source name; Col. (2): Photon Index (fixed), taken from the absorbed power-law model fits; Col. (3): The BALQSO 
subclassification; Col. (4): Ionisation parameter in erg cm s -1 ; Col. (5): Hydrogen column density of absorber in units of 10 22 cnT 2 . 
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Table 4. Hardness ratio analysis results 
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Notes: Col.(l): Source name; Col. (2): Neutral hydrogen column density estimated from hardness ratio. Errors correspond to lcr confidence 
intervals and estimated from the count-rate errors. 



